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We propose a method to create a vortex in BEC utilizing the spin degree 
of freedom. We consider the optical plug at the center of the vortex, which 
makes the vortex- creation process stable. We also investigate the instability 
of the halfway state to complete vortex state without the optical plug. 
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1. INTRODUCTION 

The vortex formation in Bose -Einstein condensate of atom gases has 
been achieved in several groups BB'EoBDQ in variety of methods. The first 
method is called phase imprinting § and an another method rotate the po- 
tential which confine the condensate, el We proposed another method in Ref. 
||, which utilized the turning magnetic field. This process is made stable 
using pinning potential with dipole force at the center of the condensate. 
This makes the condensate ring shaped and the refined properties of vor- 
tices around the core is hidden. So it is more desirable to form the vortex 
without the pinning potential. In this paper, we compare the stability of con- 
densate with and without the pinning potential. We analyze the excitation 
spectrum to obtain the insight to the stability of the process. 

The next section introduces Bogoliubov equations. Section 3 is devoted 
to the explanation of method to create the vortex state. In Section 4, we 
discuss the instability of vortex creation process. Section 5 is devoted to 
summary. 
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2. BOGOLIUBOV EQUATIONS 

In order to treat an interacting Boson system, we use the formulation 
based on Bogoliubov theory, which is extended to treat the spin degrees of 
freedom. The wavefunction of condensate is written with a set of 3 wave- 
functions (fii(r)(i = 0, ±1). These are given by the Gross-Pitaevskii (GP) 
equation 



''^^ +V(r)+g n J2\<t>k\ 2 )S ij -B ij 
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where [i is the chemical potential. The coefficient g n and g s are interaction 
parameters. The numbers Bij are matrix elements of i3-matrix which shows 
the magnetic field. The form of this matrix is given below. The subscripts are 
a = (x,y,z) and i,j,k,l = (0, ±1). V(r) is the spin-independent potential 
which is given in Eqs. (U) and ©. 

The excitation levels e q and corresponding wavefunctions u q ,v q are given 
by Bogoliubov equations 

J2{AijU q (r,j) - BijV q (r,j)} = e q u q (r,i), (2) 
3 

^2{B* jUq (r,j) - A* jVq {r,j)} = e q v q (r,i) (3) 

3 

where 

Aij = y — — -V + V(r)+g n ^ k \ 2 j6 tJ - Bij 
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Bij = g n (pi4>j + g s ^2^2(Fa)ik<t>k(F a )ji(t)i- (5) 
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Most of e q are positive. When the lowest e q is negative, it shows the local 
instability of condensate. 



3. PROCESS TOWARD THE VORTEX 



We consider a two dimensional system of BEC with rotational symmetry 
around the z axis. We introduce the cylindrical coordinates r = (r,8,z). 
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Suppose that an Ioffe-Pritchard field 

B = (f?j_ cos 9, —B± sin 9, B z ) (6) 

is applied to the system. The coefficients B± and B z are amplitude of 
magnetic field measured by Zeeman energy. This magnetic field is rewritten 
in matrix form as 



(7) 



The eigenvectors of B show the arrangement of spin components when the 
field is strong enough. The eigenvalues of B are and ±B where B = 
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B± + B z . The eigenvectors corresponding to ±B are 

(B ± B z )e id \ 

±V2B ± . (8) 
(B T B z )e- ie j 

The winding number (coefficient of 9 in exponential) becomes 1,0, and — 1 
for 4>i,(f>o and </>_i respectively. When the coefficients B± and B z vary as 

B z = B z0 cos[ir(l-t/T)] (9) 
5 ± = B' ± rsm[ir(l-t/T)] (10) 

where i (0 < t < T) is time, the spin of condensate shifts following Eq. (|8|). 
If the process begins with condensate of -1 component without vortex, the 
initial arrangement of condensate is Eq. (||) multiplied by e tS . The coefficient 
of 9 in exponential of 4>-i is zero and the condensate takes the form 

/ (B + B z )e 2ie \ 

V2B L e w . (11) 
V (B-B z ) J 

This process ends with condensate of +1 component with winding 2 at t = T. 
This means that the vortex is nucleated. 

The explanation above treat only the magnetic field. We can take into 
account other terms with GP equation Eq. (jl]). This does not change the 
arrangement of winding numbers as long as the condition of cylindrical sym- 
metry is kept. 

When the condensate follows Eq. (0), the winding numbers of excita- 
tion levels are written as qe + 2, qg + 1, qg for u+i, Uq, U-i and qg — 2,qg — 1, qg 
for v + \,vo,V-i respectively. When the pinning potential does not exist, the 
modes with qg = ±1 and ± 2 can localize at the core (and have negative 
energy) because they have winding number in Uj or Vi. 
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4. ANALYSIS OF STABILITY 

The stability of the condensate can be evaluated with sign of the lowest 
excitation level, which is given by Bogoliubov equations Eqs. (0) and (|3|). We 
decide the distribution of condensate with GP equation Eq. (Q) at each time 
t, and the excitation levels with Eqs. (|2|) and (||). The whole condensate is 
confined with spin-independent harmonic potential, and the magnetic field 
which is described in Eqs. @ and ( JlO| ) handle the spin of condensate. The 
parameters are m = 1.44 x 1(T 25 kg, a = 5.5 x 1(T 9 m, g n = ^ 2a "+ a2 , 

g s = ^p " 2 ^" - The scattering length 02 is 0.75ao or 1.375ao, which means 
g s = —0.1g n and g s = +0.1g n respectively. 

The case when the pinning potential at the center of the system is 
introduced is discussed in next subsection, and the case without the pinning 



potential is discussed in subsection 4.2. The time T is supposed to be enough 



longer than the Larmor frequency, namely the adiabatic condition is satisfied. 
When the time T is too short, the spin of atoms will not have enough time 
to follow the perpendicular magnetic field. Whether the vortex nucleates or 
not is uncertain in this case. 



4.1. WITH PINNING POTENTIAL 

The whole condensate is confined with spin-independent harmonic po- 
tential radially, and the pinning potential at the z-axis of the system is 
also given. This stabilizes the vortex creation process. These two spin- 
independent potentials are written as 



m{2irv) 2 2 , tt f ~ r<1 
2^ 



V(r)='-^^r 2 + Ue, P [- T (12) 





where v = 200Hz, U = 2 x 10 3 J and ro = lfim. 

The magnetic field which is given by Eq. @ and @ are applied. The 
wavefunctions of the condensate and the excitation spectrum are calculated 
at each t, at each pairs of magnetic field coefficient B r and B' ± . Fig. |l|(a) 
shows the density distribution of the condensate at several t. Fig. [l|(b) shows 
the lowest excitation levels at qe = ±1,±2. The change of total magneti- 
zation M along the z-axis accompanies the change of t, and the horizontal 
axis of Fig. |(b) shows M/N, which is J(\(f>i\ 2 - |</>_i| 2 )(fr/ Ei / \<t>i?dv. It 
is found that the lowest e never becomes negative. Because the negative 
excitation level does not appear, the condensate stays locally stable. 
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Fig. 1. In (a) and (c), the solid line shows the density profile of condensate 
| <j>i | 2 and the various dashed lines show the density of three spin com- 
ponents |</>i| 2 (i = 0, ±1). In (b), (d) and (e), the lowest eigenenergy e at 
winding number index qg = ±1, ±2 are plotted. The unit of e is 200/iJ, 
which is the trap unit when the pinning potential is not exist. The horizon- 
tal axis show magnetization along z-axis which varies from -1 to +1 during 
the process, (a) and (b) show shift of density distribution of the conden- 
sate and the lowest eigenenergies when the pinning potential exists. The 
eigenenergies never become negative, (c), (d) and (e) show those when the 
pinning potential does not exist. Some of eigenenergies become negative and 
even complex as the magnetization increases. In (a), (b), (c) and (d) the 
interaction parameter g s = —0.1, which means ferromagnetic interaction. a 2 
is 0.75ao. In (e), the interaction parameter g s = +0.1, which means anti- 
ferromagnetic interaction and a2 is 1.375ao- The density distributions (a) 
and (b) hardly change when the sign of interaction parameter g s is turned 
to antiferromagnetic. 
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4.2. WITHOUT PINNING POTENTIAL 

In this subsection, we consider the case when the optical pinning does 
not exist. The spin-independent trapping potential written in Eqs. (|l]) and 
©is 

V(r) = \ ) r\ (13) 

and the time-dependent magnetic field which is given in Eqs. © and © is 
also applied. The halfway states are obtained by the static form of Gross- 
Pitaevskii equation. Figure |l](c) shows some of the density distributions of 
condensate. Figure [j](d) and ffl(e) show the excitation levels obtained through 
Bogoliubov equations Eqs. (||) and (||). The filled dot shows that there is 
a pair of eigenstate with complex eigenenergy. The real parts of them are 
plotted. The complex eigenenergy appears where the qg = 2 mode and the 
q@ = — 2 mode crossed, which result is consistent with Ref. ||. This fact 
means the spiraling out or the splitting of the vortex. 



5. CONCLUSION 

In this paper, we analyzed the stability of vortex creation process. The 
distributions of condensate and excitation levels are derived using Bogoli- 
ubov equations. When the pinning potential is absent, negative and complex 
eigenenergies are found and the vortex forming become unstable halfway. 
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